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RANDOMLY FLUCTUATING FLOW IN A CHANNEL DUE TO RANDOMLY 


FLUCTUATING PRESSURE GRADIENTS 
by Morris Perlmutter 
Lewis Research Center 

SUMMARY 

A randomly fluctuating pressure gradient of a stationary Gaussian Markovian form 
will cause a randomly fluctuating velocity to be superimposed on the steady incompres- 
sible flow in a channel. Correlations, spectra, and frequency response functions for the 
random functions are given. Random pressure signals were generated using Fourier 
series expansion with coefficients randomly picked from distributions whose parameters 
were obtained from the spectra of the pressure signal. The random velocity signals 
were then obtained from the pressure signal by use of the frequency response function 
calculated from the equation of motion. The increased power loss due to the fluctuations 
is given, and the random pressure and velocity signals are compared for amplitude, fre- 
quency, and time lag. 


INTRODUCTION 

In many practical situations, flows in channels will have random pressure gradient 
fluctuations. It is of interest to understand the effect of the randomly fluctuating pres- 
sure gradient on the characteristics of the flow. 

Previously, there have been several analyses on the effect of sinusoidally varying 
pressure gradients on channel flows (refs. 1 and 2). Only recently, however, has there 
been an increasing interest in systems with randomly fluctuating characteristics. Some 
recent work on this subject was on heat conduction in solids with randomly varying ex- 
ternal temperatures or internal heat generation (refs. 3 and 4). An analysis of the effect 
of randomly fluctuating fluid velocity in a heated channel on the wall heat transfer is 
given in reference 5. 

The purpose of the present analysis is to study the effect of imposing upon a steady 
pressure gradient and steady incompressible flow in a channel a randomly fluctuating 
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pressure gradient component. Since the fluid is considered incompressible, the pres- 
sure gradient is instantaneously imposed throughout the channel. This random pressure 
gradient fluctuation will cause a random fluctuating component of velocity to be imposed 
on the steady velocity in the channel. The analysis is restricted to channel regions suf- 
ficiently downstream from the channel entrance so that entrance effects can be neglected. 
The analysis gives the mean square value of the velocity fluctuation across the channel 
as well as the increased power loss in pumping the fluid due to the fluctuations in the 
velocity. 

When a random sampling or Monte Carlo approach is used, a specific pressure 
gradient signal can be generated from a given power spectrum. This calculation can be 
rapidly carried out by use of the fast Fourier transform method of calculation recently 
developed by Cooley and Tukey and discussed in references 6 to 8. The randomly vary- 
ing pressure gradient signal is used to find the corresponding fluctuating velocity signal 
at various positions across the channel. This can be calculated through use of the mo- 
mentum equation. These results can be useful, for instance, in numerical convective 
heat-transfer calculations in which numerical values for the randomly fluctuating veloc- 
ity are required as input to calculate wall heat transfer. By generating a random signal 
with a given spectrum by the use of model sampling, these signals can be used as input 
for numerical solution of nonlinear equations. The resulting numerical output signal and 
spectra can then be found numerically. These output results can be analyzed numerically 
for such statistical characteristics as rate of zero point crossing. Solutions to nonlinear 
problems with random input are very difficult to obtain by the more usual analytical pro- 
cedures. 


SYMBOLS 


a , b complex coefficients 

n’ n 

d spacing between parallel plates 

f frequency 

f Q frequency increment, 1/T 

f maximum frequency 

H frequency response function 

| H | gain factor 

h weighting function or impulse response function 

i unit imaginary number 

K k77 


2 



N number of samples 

P pressure gradient 

R X y correlation (x, y) 

S mean square power spectral density 

T time 

Tp period of record 

2 

t dimensionless time, pT/d 

At time between sample 

U fluid velocity 

2 

u steady dimensionless velocity, -up U/d P 

s s s 

2 

v transient dimensionless velocity, -up U./d P 

L S 

Y coordinate across channel 

Y r 4 sin Ky/K 

y dimensionless coordinate across channel, Y/d 

y dimensionless pressure gradient, P^/P g 

A measure of rate of fluctuation of pressure signal 

X dimensionless measure of rate of fluctuation of pressure signal, Ad /u 
v kinematic viscosity 

p density 

2 

variance measure of pressure fluctuations 
r dimensionless time difference, tg - t^ 

cp phase factor 

<j 0 angular frequency, 2nf 

( ) ensemble average 

* complex conjugate 

Subscripts: 

F fundamental frequency or lowest frequency, 1/T 

P 

H related to frequency response function 

I imaginary 



k finite spatial sine transformed variable 
R real 

s steady component 
t fluctuating component 
y related to pressure signal 

u> Fourier time transformed variable 


ANALYSIS 

For viscous incompressible flow between parallel plates with constant properties in 
the fully developed flow region (see fig. 1) the momentum equation can be written as 
(ref. 2) 


— = - I P( T) + „ (1) 

9T P 

We can let the velocity and the pressure gradient consist of a steady part and a nonsteady 
part 


u = u s + u t 


P = 


P„ + P* 


Then equation (1) can be rewritten as 




8U t P, 

+ V 

0T P 


(2a) 

(2b) 


(3a) 


(3b) 
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Normalizing the variables as follows 


vP U fpU, 

s v = ' 

d2p s d2p s 



we obtain for equations (3a) and (3b) 



u„ 



y 


2 


(4a) 


av a 2 v 

= y + 

9t ay 2 


(4b) 


There is both a particular and a homogeneous solution to equation (4b). Since we are 
interested in times after the initial transients have decayed, we need not be concerned 
about the homogeneous solution and we can use Fourier transforms to obtain the solution 
to equation (4b). The Fourier transform (FT) pair can be written as 


v 


co 



(5a) 


v 



where 


go = 2 7Tf 


(5b) 


Taking the FT of equation (4b) gives 


3 2 v 

= y co + 2 

ay 2 


CO 


(6) 


Since the velocity must be zero at the walls of the channel, we can solve equation (6) by 
using finite sine transforms 


v 


CO 


k v “ k 


sin k^-y 


(7a) 
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(7b) 



sin kffy dy 


Then equation (6) becomes 


iwv w k = ^k- 



( 8 ) 


where 


y o)k 



4 y 

sin k7ry dy = — — 
K 


k=l, 3, . . . 


and 


K = ktf 


By solving for and taking the inverse finite Fourier transform, we obtain 

v = y H (9) 

<x> ' u> ' ' 

where 

oo 

H = V' — = H r + iHj = | H | e 1(PU (10) 

k=l^ 5 r2 + iw 

where Y K = 4 sin Ky/K. The H is sometimes called the frequency response function 
(ref. 9, p. 42), |h| the gain factor, and the phase factor, where 

|h| 2 = H 2 +H 2 (11a) 


cp H 


- tan 



(lib) 


and 


6 



(12a) 


H R = 


y y k k 

io c K 4 + w‘ 


k=l, 3, 5 


Hj = -cu 


OO 

z 

k=l, 3 


K 


K 4 + a ) 2 


(12b) 


The gain factor is plotted in figures 2 and 3; the phase factor is plotted in figure 2 
and is discussed more fully later where used. Notice that equation (9) can be written as 




v u = l v col e 


|H||r w l e 


-i i<P B +<P y ) 


(13) 


This shows that the amplitude of the Fourier transform of v is given by the product of 
the absolute values of H and y with a phase shift given by the phase angle + <Py- 
We can solve equation (9) for v by taking the inverse transform. Using the convolution 
theorem (ref. 10, p. 26) on equation (9) gives 


v 


= e iwt v H df = f + ™ h(0, y)y(t - 9)A9 

•s on ^ *S-00 


(14) 


where h is the weighting function given by 


r+°° 

h(0) = / 

‘'.on 


He iw0 di = J 


E 

k=l, 3, 5 


Y K 6 


-K 2 (0) 


for 9 > 0 


(15) 


j^O for 9 < 0 

The h and H are a Fourier transform pair. We can now write the equation for the 
fluctuating velocity as 


v(.) ./* 
‘ / 6=- 


h(t - 6, y)y(6)d6 


(16) 
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or 


/OO 

w _o h ( 0 > y)^ - 0 ) d0 


( 17 ) 


The cross correlation of the pressure and velocity fluctuations can be obtained as 
follows. From equation (17) we can write 

/ oO 

- e)de (18) 

Since the process is stationary, by taking ensemble averages, we can write equation (18) 
as 

/ OO 

h(0, y)R yy (r - 0)d0 (19) 

where tg - ti = t- Since the cross spectrum is the Fourier transform of the cross cor- 
relation, by using the convolution theorem as before, we can write 

S^M = H(w)S yy (co) (20) 


We can find the autocorrelation as 


/ OO 

h(0,y)y(t 1 - 0)v(t 2 )d0 

Ensemble averaging as before gives the velocity autocorrelation 

RwW = h(0, y)R yv (r + 9)d9 

We can write the power spectrum as 


( 21 ) 


( 22 ) 


8 



/'+ co 

s w (w) = L 


-lOJT 


/’CJO 

/ h(0, ylR^fr + 0)d0 

J d=0 7 


dr 


= f h(0, y)e ic,;0 dd f + R ( T + 0)e' ia,(T+0) d( T + B) (23) 

^0=0 ^°o ^ 


This gives 


S vv (co) = H+MS^M = H*(co)H(co)S ry (a;) 


(24) 


S vy (co) = |H(co)| 2 S yy (co) (25) 

This shows the velocity spectrum is equal to the pressure spectrum times the 
square of the gain factor. 


Pressure Gradient Correlations and Power Spectrum 


A common form of random fluctuations encountered in physical systems is called a 
stationary Gaussian Markoff random process. Markoff processes are random processes 
whose relation to the past does not extend beyond the immediately preceding observa- 
tion. As is shown in reference 11, page 215, and reference 12, the autocorrelation 
function for this random process is given by 


R = ct 2 e A l 

77 y 


= a 2 e- X 

r 


(26) 


where 


< 4 = v o) s <y2) 


(27) 


The rate of fluctuation A and the time span T have been nondimensionalized to be 

j2 A 


A d 


> 


(28) 


T = 


_ Tv 
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The Fourier transform of gives the power spectrum of the pressure fluctua- 

tions as a simple Lorenzian form: 



2 Xa 


y 


.2 2 

X + w 


(29) 


2 

The a is a dimensionless measure of the amplitude of the random fluctuations 

2 y 2 

(y ) . The larger the value of cr^, the larger the amplitude of the pressure fluctuations. 
The A. is a dimensionless measure of the average number of fluctuations per unit time. 
It can also be considered as the inverse of the dimensionless characteristic decay time 
of the autocorrelation. The larger the value of X, the greater the fluctuation rate of 
the signal. The power of the pressure signal is given by the integral of the power spec- 
trum over all frequencies. This can be written as 


r (o) = / S df = ct 2 = (y 2 ) (30) 

yy J-oo yy y 

9 

The normalized power spectrum for the pressure gradient fluctuations S ^/(y ) is 
plotted in figure 3 for different values of the rate of fluctuation parameter X. It can be 
seen that for smaller values of X much more of the spectral power is located in the 
lower frequencies since the frequency band for X of 0. 1 is much narrower than in the 
case for X = 100. 


Fluctuating Velocity Mean, Correlations, and Power Spectrums 

The pressure-velocity cross correlation for the present case can be written using 
equations (15), (19), and (26) as 
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This can be integrated to give 


r 




V t) = < 


•si 

k=l, 3, 5 

CO 

<z 


L K 


e-X T _ __2X_ e -K* T 


K - X L X + 


for t > 0 


> (31b) 


Xt 


for r < 0 


v 


k=l, 3, 5 


K 2 + X 


Also the cross spectrum can be written using equations (10), (20), and (29) as 

.2 


V w) 


I 

Kk=l, 3, 5 


y k \ 2x % 


„2 1,2 2 

K + ico / X + Cl) 


(32) 


The result in equation (32) can also be obtained by taking the FT of equation (31). 

In a similar manner we can calculate the velocity autocorrelation from equations 
(15), (22), and (31) as 



(33) 


This can be integrated to give 


R (t) = o 
w v ' ' y 


oo OO 



k=l, 3 £=1,3 


y k y l 

e~ Xr , 

( 2X ^ 

T 2 1 

e- L t 

l 2 -x 

_X + K 2 

vx + l 2 ; 

'k 2 + l 2 


(34) 


Similarly, we can write the velocity fluctuation power spectrum using equations (12), 
(25), and (29): 


SwM = |H(co) | 2 S yy (u>) = (H 2 + 


(35) 
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From equation (16) we can see that the ensemble average of the fluctuating velocity 
(v) at any given y is zero since the ensemble average of the fluctuating pressure gradi- 
ent (y) is zero. 

o 

The mean square velocity (v ) can be obtained from equation (34) since 

R w(°’ y ) = < y2 > ( 36 ) 


These results are shown in figure 4. Notice that the mean square velocity is equal to 
the total power of the velocity spectrum 


/ +°° /'+°° 9 

S vv df =J oo iHM^Mdf 


(37) 


The power or the mean square velocity is given by the product of the square of the gain 
factor times the pressure spectrum integrated overall frequency. From figure 2 we 
can see that the magnitude of the gain factor is smaller near the wall (y = 0. 1) compared 
to the centerline (y = 0. 5), so that the mean square velocity is smaller near the wall as 
compared to the centerline. The mean square velocity goes to zero at the wall in fig- 
ure 4, as would be expected from the zero wall velocity boundary condition. The gain 
factor in figure 2 is greatest at the lower frequency and has a high-frequency cutoff. 

Since the spectrum for the pressure gradient (fig. 3), is also greater at the lower fre- 
quencies for low values of X as compared to the high X case, the mean square velocity 
will be larger for the lower values of X, as can be seen in figure 4. 

The fluctuating velocity power spectrum given by equation (35) normalized by the 
total power of the spectrum at a given position in the channel S yy (y, f)/R yy (y, o) is also 
shown in figure 3. These curves are discussed more fully later. 


Power Dissipation Due to Velocity Fluctuations 

It is desirable to know how much additional power must be supplied to the flow sys- 
tem to maintain the same average flow rate for the fluctuating flow as compared to the 
steady flow case. As shown in references 1 and 2, the external rate of work done by the 
pressure force w g is equal to the sum of the rate of increase in kinetic energy K 0 and 
the rate of dissipation of energy due to internal friction w^: 

dK 

w = — - + w f (38) 

e dT 1 


12 



This can be found also by multiplying the momentum equation (1) by velocity and inte- 
grating across the channel to obtain 


-UP = 



dy 


( 39 ) 


where U is the average across the channel of some quantity. 

If we then take the ensemble average, the rate of change of the kinetic energy term 
will be zero because the process is stationary. We then have the external rate of work 
done by the pressure force w g equal to the rate of dissipation of energy due to external 
friction w^ 


-(UP) 



(40) 


To find the rate of dissipation of energy due to internal friction we can then write 


W { = -( UP) = "UgPg - ( U t P(> - w tj , - < U t P t > (41) 


If we let w f be the rate of dissipation due to steady flow, the ratio of increased dis- 

I, fa 

sipation due to the unsteady flow over the steady flow dissipation can be written as 


f - w f <U,P.) _ _ 

^ = 12( vy) = 12R (o) (42) 

w f, s u s p s y 


Integrating equation (3 lb) over y and substituting into equation (42) give 


Wf ~ Wf » s = 96a 2 V * (43) 

Wf > s 7 feW 5 r2 ^ 2 + 

A plot of equation (43) is given in figure 5 as a function of the fluctuating rate param- 
eter As X becomes very large, that is, fluctuations of the pressure gradient be- 
come very rapid, the frictional power loss reduces to the steady power loss. The fric- 
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tional power loss increases with smaller values of X or slower fluctuations of the 
pressure gradient. 


Generation of Random Signal by Model Sampling 

It would be helpful in certain cases to be able to generate a random signal analyti- 
cally on the computer consistent with the power spectrum of the signal. This would be 
useful, for instance, in obtaining numerical solutions or various other statistical char- 
acteristics of linear or nonlinear equations of random systems that may be difficult to 
obtain by more usual analytical procedures. The random signal can be obtained by 
writing the signal in terms of a Fourier series expansion. The coefficients of the ex- 
pansion can be related to a random distribution whose parameters can be obtained from 
the power spectrum. We can then use random sampling to choose the coefficients from 
a distribution with the appropriate parameter and thus obtain the coefficients of the 
Fourier expansion. We can then make use of the fast Fourier transform method (refs. 

6 to 8) to solve for the signal given the coefficients. The details are given in the next 
section. 


Relation Between Power Spectrum and Fourier Transform 

First we must relate the power spectrum to the Fourier transform of the random 
signal. This can be done as shown in reference 12, page 132. We can define a function 
Z Yy {T) given by 


i r+ T/2 

# vv (t) = lim - J y(t)y(t + r)dt 
yy t—°° t -t /2 

Using Parsevel's formula and the time shift theorem (ref. 10, pp. 14 and 27), we can 
write 


9t- At) = lim I f 
* > T— oo T — °o 


+00 


* „iU>T jr 
y y e dl 
'or u> 


(45) 


Taking the ensemble average of equation (44) gives 

•T/2 


{ 0 yy > 2 V T)dl = V T) 
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Taking the ensemble average of equation (45) gives 


R (t) = lim 
yy ' T _oo T 






( 46 ) 


Since the power spectrum and correlation form a Fourier transform pair, we can write 




(47) 


This, then, relates the power spectrum to the mean square of the Fourier transform 
of the random signal. 

Discrete Fo urier transforms . - When a waveform is to be analyzed on a digital com- 
puter, it is the discrete Fourier transform that must be used. This can be derived from 
the usual Fourier transform as shown in reference 11, page 56, as follows. 

If y(t) is periodic in T , we can write the Fourier expansion 

Jr 


y(t) 




in27rt/T 


n=- c ° 


(48) 


where 


1 /V 2 

a n=“ / 

p«AT /2 


-i27rnt/T 

y(t)e p dt 


(49) 


We can see from equation (48) that 


X 11 “„ (T P a n> = 

P 


(50) 


Similarly, if y is periodic in f we can write 

Ct) p 




= Tj 6 ‘ 


27rnf/f 


n=-°o 


(51) 
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where 



We can see that, similarly, 

lim ( f n d n ) = y( 4 ) (53) 

f —OO V 11 

p 


If we take the time between sample values of y k as At, we can take N samples 

so that N At = T . Then, assuming that At = l/f we can write equation (48) as 
P P 


N/2 

y k (k At) = £ a e 

K n=^N/2 n 


i27rkn/N 


(54) 


We can also write equation (51), assuming that N Af = f as 


N/2 


a n = — E, y k e- i2 ' k "/ N 

n N k=TN/2 K 


(55) 


Since a n can be related to the Fourier transform y^ by equation (50), equations (54) 
and (55) form a discrete Fourier transform pair that is an approximation of the Fourier 
integral transform pair. It can be shown (ref. 6) that in fact equations (54) and (55) are 
a discrete transform pair with reversible mapping in its own right. 

Fast Fourier transform. - A numerical method for evaluating either equation (54) 
or (55) in a very efficient manner on digital computers has been developed called the 
fast Fourier transform (refs. 6 to 8). One implementation of the fast Fourier transform 
is available as an IBM- scientific subroutine package called HARM/DHARM. Many 
others are also available. 

Use of the fast Fourier transform allows the spectrum to be found by taking N 

samples of a given signal y(t) of length T p in time increments At such that 

N At = T . Using the fast Fourier transform method, we can readily solve for the real 

( a n p) and imaginary (a n j) Fourier expansion components in equation (55). Repeating 
5 ’ k 

this process for k signals, we can obtain the averaged quantities (a^ j) -m< i), 

k i=l 1 
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o 

and similarly (a* p) . The spectrum of y(t) can then be obtained from equations (47) 
and (50) as 



w=n27r/T 

P 





(56) 


Gaussian random process . - Equation (54) can be rewritten by letting 

= a^ D + ia„ T to obtain 
n n, R n, I 


N/2 

^k = z 

n=-N/2 


(a n,P + ia n,I 


Jfcos^+isin^A 
N N / 



a„ T sin 
n, I 


in 2j[nk\ 

N J 


+ i( a„ tj sin ^ 7rn ^ + a T cos 
i n, k N n, i 


27rnk \ 
N /. 


(57) 


Since y^ is a real function, the imaginary term in the right side of equation (57) must 

be zero. If a^ D is an even function around n = 0, a „ 0 = a„ t-,, and if a„ T is an 
n, R ’ -n, R n, R’ n, I 

odd function, a T = -a T . Then we can see that 
-n, I n, I 


N/2 

I 

n=-N/2 


a n, R sin 


2?rnk 

N 


-1 

z 

n=-N/2 


a sln2lSi + 

n ) It JyJ 


N/2 

z* 

n=l 


n, R 


sin*E!* = 0 
N 


(58) 


Similarly, 



n=-N/2 


a t cos2E* = 0 
n ’ 1 N 


so that we now have 
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XT /O 


(59) 


We can assume a Gaussian random process (ref. 12, p. 155), that is, a process where 
the distribution of y at different times [y(tj), y^), . . .J is a Gaussian distribution. 
This process is commonly used as an idealization of many natural phenomena associated 
with the superposition of a large number of many small effects. Since a n can be seen 
as a summation of Gaussian variables from equation (55), a n must also be a Gaussian 
distribution (ref. 12, p. 78). It is shown in reference 12, page 160, that for a Gaussian 
random process the following conditions are true: 



< a n, R a k, R> ~ ^ a n, I a k, 1^ 


0 


k=n 

► 

k^n 

J 


<Vk> = 0 


Then, equation (55) can be written as 


S = 2T (af) when n^O 

^ U)=n27r/T p 


However, since a T is zero when n = 0, from equation (56) we have 
n, i 



= T 


P 


< a 0> 


(60) 


(61) 


(62a) 


(62b) 


Model sampling to genera te signal . - To generate a random signal y(t) we can use 

equation (54) and the fast Fourier transform program. The terms a „ and a T are 

n ? k n, i 

random variables given by the joint Gaussian distribution 


f < a n, R’ a n, 


f( a n, R +a n, I^ 2 ( a n) 


2^<a l) 


(63) 


Following the method given in reference 13, page 39, we can randomly choose values of 
a n R anc * a n I *kat wi ^ satisf y the distribution of equation (63) by using the following 
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equations: 


a »,R = ( 2 < a n>) R r ) 1/2 cos 2irRg 

(64a) 

a n, I = ( 2 < a n>) V2( - ln R r )1/2 6in 2 " B « 

(64b) 


where R fl and R are two different random numbers uniformly distributed between 
0 and 1. The values of R can readily be generated by the computer. The values of (a^) 
can be obtained from equations (62a) and (62b) using the appropriate power spectrum 
S . Recalling that a n ^ is an even function and a n j is an odd function, in the par- 
ticular fast Fourier program available (HARM/DHARM), we only need find a n „ and 
a n j from n = 0 to N/2. Then a^ N y 2 ) +1 R = a (N/2)-l R’ etc - Similarl y> for a n i 
values were only needed for n = 0 to n = N/2; since a (N/2)+l I e( l ua ll e d ~ a (N/2)-l I> 
etc. , a Q j was equal to zero. 


Pressure Gradient Signal Results 

, 2 1 / 2 

The normalized fluctuating pressure gradient y/(y ) is shown in figure 6 for 
different values of X. The pressure gradient power spectrum used in obtaining these 
results was assumed to be a stationary Gaussian Markoff process and is given by equa- 
tion (29) and plotted in figure 3. The results were calculated for values of N of 128 and 
XTp taken as 100. 

The curves were computer plotted by joining the output points by dotted lines. Since 
the spectrum for the higher fluctuating rate parameter X = 100 is greater at the higher 
frequencies, as can be seen in figure 3, the resulting pressure signal for large X has 
many more fluctuations and crossing of the 0 value line per unit time than is the case 
for the lower values of X. 


Generating Velocity Signal From Pressure Signal 


The method of calculating the velocity signal from the pressure signal is as follows. 
We can write, similar to equation (59), 


v 


k " 




cos 


27rnk 

N 


+ b 


n, I 


sin 2 !^ 
N / 


(65) 
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where b n ^ is an even function and b n j is an odd function. From equations (9) 
and (50) we can write 



v a>= b n T p = P I | # )( a n T p) 


or 

b n = H n a n 

(66) 

This can be written as 



b n 

= b n, R + lb n, I = ^ H n, R + lH n, I^ a n, R + ia n, P 

(67) 

The components of b n 

are then given by 



b n, R = a n, R H n, R " a n, I H n, I 

(68a) 


b J — 3. tH T-) + 3 f 

n, I n, I n, R n, R n, I 

(68b) 


At n = 0, since a^ T = H T = 0, it follows that b„ D = a^ 0 H d and b T = 0. Thus, 
o, 1 o, 1 ’ o, K o, R o, R o, I ’ 

by finding the values of a n ^ and a n j for a given signal y^, as discussed previously, 
we can use equations (68a) and (68b) to obtain the new values b n ^ and b n j, which can 
be used in equation (65) to obtain the velocity signal v using the fast Fourier transform 
program. 

We can see from equation (68) that the mean values of b R ^ and b n j are zero. 

Also, since a n and a T are normally distributed, b D and b T are also nor- 
’ n, R n, I J ’ n, R n, I 

mally distributed. The variance of b n can be seen to be 



< a n, R^ H n, R 



Since 


< a n,R> equals (a^j), 


it can be seen that 


(69) 


<<R> '<<!> '< a n>( H n,R* H n, l) = < a n> I H n 1 2 


(70) 



except at n = 0 when 




/h 2 s _ ,2 x „2 

^ b o, _ < a o, R^ H o, R 




J 


The spectrum in this case would be given as in equation (56) by 


S v>> = T p 2 < b n> = T p 2 < a n>l H | 2 n *° 


so that 


2T <b 2 ) = S (n) |H I 2 
p\ n 7 yy ' n 1 


(71) 


(72) 


(73) 


except at n = 0 when 


T <b 2 > = S (o) I H I 2 
p\ o' yy v ' 1 o 1 

This shows that b R and b T cannot only be found from a but can be picked 
iij it n, j. n n 

directly from a normal distribution with a variance of (b^) given by equation (73) as in 
the case of a n values. 

In figure 6 the normalized fluctuating velocity at the centerline of the channel 
/( 2 ^ 2 \ 

v y=0 5/v v ) J y-0 5 is P^ ot ^ ec ^ as s °lid lines. These results were calculated from 
equation (65) using the fast Fourier transform routine. The coefficients b n were found 
from equation (68), which gives the values of b R in terms of values of a n which are the 
coefficients of the random pressure fluctuations. These were randomly chosen by equa- 
tion (64). Thus, the velocity shown by the solid line on each graph is caused by the pres- 
sure fluctuation shown as the dotted line on the same graph. Notice that all the curves 
have been normalized by their own root mean square values. 

For low values of X, the shapes of the curves are very similar but there is a small 
time lag in the velocity curve compared to the pressure curve. This time lag can be 
explained by postulating a distortionless frequency response function H’ with a distor- 
tionless gain factor |H(co)| = |h| t and a distortionless phase factor cp' = -^t Q (see 
eq. (11)). 
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Then equation (9) can be written as 


icot 




H 


y, 


<J0 


(74) 


Taking the inverse transform and using the time shifting theorem (ref. 10, p. 14) gives 


V(t) = | H ' y(t + t Q ) 


(75) 


This shows, for the distortionless response function, the velocity signal lags the pres- 
sure signal by a time increment t , while the amplitude is changed by an amount |H | '. 

In comparing the distortionless frequency response function to the present response 
function we can see from figure 2 that for small angular frequencies the response func- 
tion approaches the distortionless case. We can calculate the time lag t for small 
angular frequencies by solving for cp as co — 0. This gives, from equation (11), 

(p = -^(0. 0781). Thus, the time lag between the velocity and pressure signals for 
small w is given by t = 0.0781. Thus, in figure 6 for A=1 we see the velocity sig- 
nal is close to the pressure signal with a small time lag between them. This is because 
the response function is close to distortionless for small u> and the pressure signal 
for small A is concentrated in the small u> region (see fig. 3). 

For the higher values of A we can see that the velocity signal does not contain many 
of the higher frequency fluctuations that are present in the velocity signal. This is be- 
cause the gain factor (see fig. 3) has a high-frequency cutoff which caused the velocity 
spectrum to have a similar high-frequency cutoff. Since at high X much of the pressure 
spectrum is in the higher frequencies, these higher frequencies are not transmitted to 
the velocity signal. 

In figure 7, the velocity signal at a given position y in the channel, normalized by 
its root mean square value at the same y value, is plotted at two different values of y 
for the same pressure signal. For low values of X the curves near the center of the 
channel (y = 0. 5) are very similar to that nearer the wall (y = 0. 01). This happens be- 
cause the transfer function for small u> is close to distortionless. The velocity signals 
near the wall, although of a smaller amplitude than the velocity signal near the center of 
the channel (see normalizing factor in fig. 4), both have almost identical shapes. At the 
higher values of X we see that the centerline signal fluctuates much less than the veloc- 
ity signal near the wall. This is due to the fact that the gain factors for y close to the 
wall are flatter and do not have as strong a high-frequency cutoff as the gain factors 
near the center. This causes the velocity spectrum near the wall (y = 0. 1) given in fig- 
ure 3 to be flatter and extend to higher frequencies than the velocity spectrum near the 
center of the channel. 
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At the higher values of A. = 100, we can also notice a time lag between the signals. 
This effect is due to the phase factor shown in figure 2. It can be seen that the phase 
factor curves for the signals closer to the wall have smaller slopes and so have a 
smaller delay time than the velocity near the center of the channel. 


RESULTS AND CONCLUSIONS 

The effect of superimposing a randomly fluctuating pressure gradient of a stationary 
Gaussian Markovian form with a zero mean value upon a steady pressure gradient will 
cause a randomly fluctuating velocity component to be superimposed over the steady 
flow in a channel. The mean value of the velocity fluctuation is zero. The mean square 
value of the velocity fluctuations which is a measure of the amplitude of the fluctuations 
will be highest in the center of the channel and reduce to zero at the wall. Also, the 
slower the rate of the pressure fluctuations, the greater the amplitude of the velocity 
fluctuations. 

The frictional power loss in pumping the fluid is increased by the pressure fluctua- 
tions, the slower pressure fluctuations giving larger power losses. 

The normalized velocity signal values were very similar to the normalized pres- 
sure signal but with a dimensionless time lag t of about 0. 0781 for low pressure fluc- 
tuation rates. At the higher pressure fluctuation rates the higher frequencies of the 
pressure signal did not appear in the velocity signal. 

Comparison of the velocity signals across the channel normalized by its local root 
mean square value showed that for low fluctuation rates in the pressure signal the values 
of the normalized velocity near the wall were almost identical to the normalized veloc- 
ities near the center of the channel. However, at the higher pressure fluctuation rates 
the higher frequency fluctuations present in the velocities near the wall do not appear in 
the velocities near the center of the channel, and the velocities near the center of the 
channel lag the velocities near the wall. 

The random pressure signal function was generated by a Fourier series expansion 
where the coefficients were randomly chosen by model sampling from a frequency dis- 
tribution whose parameters are given by the power spectrum of the signal. By use of 
the fast Fourier transform method of computation the function can be rapidly evaluated 
from its Fourier coefficients. The Fourier expansion of the fluctuating pressure signal 
can then be used in the momentum equation to calculate the randomly varying velocity at 
various positions from the wall. These velocity results can be useful when the details 
of a randomly fluctuating velocity distribution in a channel are needed, as for instance 
in convective heat-transfer problems. 

This method can also be useful in analyzing other linear and nonlinear differential 
equations of systems with random behavior. By generating the input random signal from 



a given spectrum, the differential equations can be solved numerically to give the output 
random signal, which then can be studied numerically to obtain its statistical behavior. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, November 24, 1970, 
129-01. 
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Figure 1. - Parallel plate channel flow model. 
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Figure 5. - Percentage increase in power dissipation ratio due to flow fluctuations. 
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(a) Dimensionless pressure gradient fluctuation rate, 1. 



Figure 7. - Fluctuating velocity signals at different positions in channel 
normalized by their local root mean square value. 
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Normalized fluctuating velocity, 




(d) Dimensionless pressure gradient fluctuation rate, 1000. 
Figure 7. - Concluded. 
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